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Phototransistors are attractive devices for applications in optical fiber 
telecommunication systems. They are used for the detection of optical signals 
and the amplification of these signals. This paper presents an investigation of 
how the technological parameters of indium phosphide (InP)/indium gallium 
arsenide (InGaAs) heterojunction bipolar phototransistor can impact its 


responsivity at two wavelengths, 1310 nm and 1550 nm. Based on the results 
of this investigation, we proposed optimized structures for the studied 
phototransistor. In this work, we used the software technology computer 
aided-design (TCAD)-Silvaco to simulate the physical and the electrical 
behavior of the different structures. The proposed optimized phototransistors 
can be used for various optoelectronic applications. 
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1. INTRODUCTION 

Semiconductor devices are indispensable for communication and information systems, they are more 
requests for various applications. III-V semiconductor materials are used in the fabrication of electronic devices 
owing to their transport properties. They have a direct band gap, and their electronic mobility is high. They have 
exceptional speed characteristics, furthermore they enable the operation of electronic devices at high frequencies [1]. 

Indium phosphide (InP) and indium gallium arsenide (InGaAs) are semiconductor materials of the II- 
V family of semiconductors. InP-InGaAs HBTs have proved to be excellent electronic devices, they are 
characterized by high electrical performances of speed and superior frequency [2], [3], and excellent current 
handling capability [4]. The heterojunction bipolar phototransistor (HPT) has almost the same epitaxial layer as 
the HBT, but the exclusive difference is in the electrode structure which allows illumination at the top of the 
device structure [5], [6]. The optical signal is absorbed in the base-collector region, and the signal is amplified 
by the intrinsic gain of the transistor [7], [8]. 

Different research works on HPT were performed for several decades [9]-[11]. HPTs are considered 
promising alternatives to devices such as metal-semiconductor-metal (MSM), the PIN and avalanche 
photodiodes [12], [13]. Furthermore, HPTs are more advantageous for highly sensitive infrared photodetection 
and imaging over the PIN and the avalanche photodiode (APD) due to their excellent characteristics of high 
gains without differences in high potentials and without excess noise [14]-[16]. In addition, the operation of 
HPTs can be done at a low bias voltage, and with a low gain sensitivity to the bias voltages [17]. There are 
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several applications where photodetectors based on InP are requested and used, such as biomedical field, optical 
communication systems, environment and also firms of chemistry for near infrared (NIR) light detection [18]. 

The objective of this work is to investigate the impact of technological parameters on the electrical 
performances of the InP/InGaAs HPT in terms of the responsivity at the wavelengths 1310 nm and 1550 nm. 
This investigation will help us to propose some optimized devices with different values of responsivities. These 
devices can be used in various optoelectronic applications. 


2. INDIUM PHOSPHIDE/INDIUM GALLIUM ARSENIDE HPT MODELLING 

In this work, we designed the InP/InGaAs HPT, which is mainly based on the HBT structure studied 
in papers [2], [19]—[21]. It is basically composed of two IIJ-V semiconductor materials made from two alloys, a 
binary alloy for InP, and InGaAs. The topology of HBT structure has been adapted to receive light. Regarding 
the modifications made in HBT structure, a base contact has been removed, and the emitter contact has been 
reduced. The type of illumination is from above. 

As detailed in Table 1 and shown in Figure 1, the epitaxial layer structure of the InP/InGaAs HPT 
comprises a cap layer, an emitter | layer, an emitter 2 layer, a spacer layer, a base layer, a collector layer, a sub- 
collector layer and a buffer layer. For the fabrication of this device, there are different techniques used for the 
growth of epitaxial layers such as metal organic chemical vapor deposition (MOCVD), and molecular jet epitaxy 
(MBE) technique [22]. 

Figure | shows the device structure of the InP/InGaAs HPT that respects the order of the epitaxial 
layers of the Table 1. It is an NPN heterojunction bipolar phototransistor; the ohmic contacts are made of the 
material gold. In the top view, the optical window as illustrated is the space between the base and the emitter. 
The illuminated surface is equal to 4x4 m?. The illumination is considered vertical to the HPT optical window. 
To design the device structure, we have carefully layered and meshed the InP/InGaAs HPT using the device 
editor DevEdit which allows the creation of the device structure. 


Table 1. Layer structure for lattice-matched InP/In9.47GQ9.53 HPT 


Layer Material Doping (cm=?) Thickness (nm) 
Cap Ino.47Gag.s3As n=1x101 135 
Emitter 1 Inp.47Gao.s3As n=1x1077 135 
Emitter 2 InP n=1x10"” 40 
Spacer Ing.47Gag53As - 5 
Base Ing 47Gags3As p=1.5x10?9 65 
Collector Ing.47Gags3As n=1x10"° 630 
Sub-collector Ing 47Gags3As n=1x1019 500 
Buffer Ing.47Gag53As - 10 
Substrate Semi-insulating InP 
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Figure 1. Schematic device structure of the InP/Ing 47Ga9.53 HPT 


2.1. Physical and numerical modelling of InP/InGaAs HPT 
ATLAS is a device simulator. It helps users predict the electrical behavior of any electronic device. It 
has been integrated into the simulation program for the physical and numerical modelling of HPT. 
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2.1.1. Physical modelling 

The physical modelling of electronic devices allows to take into consideration the physical phenomena 
and mechanisms occurring within these devices. After designing the HPT device structure, we added in the 
simulation some physical models integrated into the simulator. Among these physical models, we cite the 
carrier statistical model (BGN), the recombination model (SRH), the Selberherr's model of the impact of 
ionization IMPACT SELB), the band-to-band model (BBT.STD), the optical recombination model (OPTR), 
the parallel electric field dependence model (FLDMOB), and the fermi-dirac model (FERMI). 


2.1.2. Numerical modelling 

Semiconductor equations were solved employing the newton method. Among these equations, the carrier 
continuity equations, the Poisson’s equation, the drift and the diffusion equations for both the electrons and the holes. 
These equations allow calculating currents, potentials and charge carriers, more details can be found elsewhere [21]. 


3. RESULTS AND DISCUSSION 

This study is carried out at room temperature T=300 K, for the two wavelengths 1310 nm and 1550 
nm. The photodiode mode is taken for at V..=0 V and V,,=0 V. The phototransistor mode is considered for the 
bias conditions at V..=1.6 V and V,,=0.86 V. 


3.1. HPT responsivity 
3.1.1. Responsivity definition 

The responsivity R (A) is defined as the ratio of photo-generated current to the input optical power and 
is expressed as (1) [5], [23]: 


Alc(A) 
he®g(1-Rf) 


RA) = (1) 


where / is planks constant, c is the velocity of light, ®p is total incident flux density and Ry is the Fresnel 
reflection coefficient. The flux-dependent collector current (/-) is the sum of individual current contributions 
in the base, sub-collector and B-C depletion regions [23]. 

The phototransistor is working in two modes, photodiode mode and phototransistor mode; we are 
more interested in the phototransistor mode. The figures below present the simulation results of the studied 
InP/InGaAs HPT. We simulated the collector current as a function of light intensity, after that we extracted the 
responsivity in phototransistor mode Rypr. For the applied light intensity, it is from 0 to 100 W/cm?. All the 
results for responsivity are taken for the phototransistor mode, more precisely for V%..=1.6 V and V,.=0.86 V. 
We investigated the impact of technological parameters on the responsivity of the studied HPT. 


3.1.2. Responsivity of the studied HPT at 1=1310 nm and 2=1550 nm 

Figure 2 and Figure 3 present the evolution of the collector current (A) as a function of the light intensity 
(W/cm?) respectively at A=1310 nm and A=1550 nm. The light intensity is from 0 to 100 W/cm?. It is clearly 
noticed a linear evolution, the collector current increases gradually with the increase of the applied light intensity. 
The responsivity Rypr is approximately equal to 3.2 A/W at A=1310 nm, and it is around 2.7 A/W at A=1550 nm. 
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We have investigated the impact of technological parameters of the studied phototransistor on its 
responsivity in the phototransistor mode. This investigation was performed at room temperature T=300 K 
respectively for the wavelength A=1310 nm and A=1550 nm. The evaluated technological parameters are the 
base width, the base doping, and the collector doping. 
A. Impact of technological parameters on responsivity in phototransistor mode at A=1310 nm 
1) Impact of base width 
Based on Figure 4, the collector current (A) increases with the application of light intensity (W/cm?). 
We noticed that the collector current increases with the reduction of the base width W,. Figure 5 shows the 
responsivity Rypr(A/W) of the HPT that varies according to the base width W, for the wavelength A=1310 nm. 
We observed that the responsivity had an increasing trend with the decrease in the base width. Therefore, the 
decreased values of base width allowed obtaining high values of the responsivity of the HPT, the intrinsic 
current gain of a transistor increases when we reduce the base width [20]. The base width has the same 
important impact for the intrinsic current gain and the responsivity. 


0.00084 — 5 
~ 0.00082 — af 
< =| Ss 4 
> 0.0008 = 
< a <q 3.5 
0) he = 
= 0.00078 = > 3 
O 0.00076 — 2 25 
i 4 oa 
£ 0.00074 = 5 2 
Fi 4 o 15 
2 0.00072 = 3° 
Fs} 5 ec 1 
0.0007 
1 0.5 
0.00068 ——, T t t 7 7 0 


0 20 40 60 80 
Light intensity (W/cm?) 


Figure 4. Collector current as a function of light 


intensity for different base width W, at A=1310 nm 
and V.,=1.6 V 


2) Impact of base doping 
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Figure 5. Responsivity as a function of base width 
W, at A=1310 nm and Y%,=1.6 V 


Figure 6 presents the evolution of the collector current (A) as a function of the light intensity (W/cm?) 
for different base doping N, at A=1310 nm. According to Figure 7, we find that when the base doping Ny, is 
reduced, the responsivity Rypr(A/W) increases gradually and reaches a large value for N, = 1x1018cm73. The 
base doping N, has a significant impact on the responsivity of the studied phototransistor. 
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Figure 6. Collector current as a function of light 
intensity for different base doping N, at A=1310 nm 
and V..=1.6 V 
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Figure 7. Responsivity as a function of base doping 
N, at A=1310 nm and V.=1.6 V 
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3) Impact of collector doping 

According to Figure 8, the collector current J, increases with the increasing of applied light intensity 
(W/cm*) at A=1310 nm. The investigated values of collector doping are between 1x10** cm™? and 
1x10*” cm~3. Based on Figure 9, we noticed from the examined values that the responsivity Rypr(A/W) 
increases slowly and its value is constant for some values of N,. Collector doping has a small impact on HPT 
responsivity compared to other technological parameters. 
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Figure 8. Collector current as a function of light Figure 9. Responsivity as a function of collector 
intensity for different collector doping N, at doping N, at A=1310 nm and V,,=1.6 V 


4=1310 nm and V,,=1.6 V 


B. Impact of technological parameters on responsivity in phototransistor mode at A=1550 nm 
1) Impact of base width 

Figure 10 presents the curve of collector current [¢ as a function of the light intensity for different 
values of base width W; between 45 nm and 85 nm. According to Figure 11, it is observed that the responsivity 
Rypr(A/W) decreases with the increasing of the base width W, at the wavelength 1550 nm. The base width as 
a technological parameter has an impact on the responsivity of this studied photodetector. 
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Figure 10. Collector current as a function of light Figure 11. Responsivity as a function of base width 
intensity for different base width W, at A=1550 nm W,, at A=1550 nm and V,,=1.6 V 
and V,.=1.6 V 


2) Impact of base doping 

Figure 12 presents the evolution of the collector current as a function of the light intensity (W/cm?) 
for different base dopings N, at A=1550 nm. Based on Figure 13, we find that the responsivity Rypr(A/W) 
increases when the base doping N, is reduced. The base doping N, corresponding to 1x101® cm~? gives the 
highest responsivity among the investigated base dopings ranging from 1x1018 cm=3 up to 1x10?° cm73. 
Therefore, the base doping N, has a considerable impact on the responsivity of the phototransistor. 
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Figure 12. Collector current as a function of light Figure 13. Responsivity as a function of base doping 
intensity for different base doping N, at A=1550 N, at A=1550 nm and V,,=1.6 V 
nm and V,,=1.6 V 


3) Impact of collector doping 

According to Figure 14, the collector current J¢ increases with the increasing of the applied light 
intensity (W/cm?) at 4=1550 nm. The investigated values of collector doping range from 1x101* cm-? to 
1x1017 cm73. Based on Figure 15, we noticed that responsivity Rypr as a function of collector doping N, 
varies in a very light way. The responsivity increases slightly. 
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Figure 14. Collector current as a function of light Figure 15. Responsivity as a function of collector 
intensity for different collector doping N, at doping N, at A=1550 nm and V,,=1.6 V 


4=1550 nm and V,,=1.6 V 


3.1.3. Proposition of optimized structures 

Based on the study of the impact of technological parameters on the responsivity of the phototransistor 
HPT InP/InGaAs. We proposed four new optimized topologies as shown in Table 2. These topologies can be 
used as photodetectors in optoelectronic applications. The Figures 16-23 present the simulation results of 
responsivity for the new optimized topologies of the InP/InGaAs HPT for A=1310 nm and A=1550 nm. 


Table 2. Proposed optimized topologies HPTs 
Studied Optimized Optimized Optimized Optimized 


Changed parameters 


structure structure A structure B structure C structure D 
Base width W, (nm) 65 50 50 50 50 
Base doping N,(cm~*) 1.5x107° 1x10?° 5x1078 2x10%8 1x1078 
Collector doping N,(cm~?) 1x107° 1x10?” 1x10?” 1x10?” 1x10?” 


A. Responsivity of the optimized structure A 

Figures 16 and 17 present the evolution of the collector current J; as a function of the light intensity in 
phototransistor mode at V..=1.6 V respectively for A=1310 nm and A=1550 nm. The collector current [¢ increases 
with the increasing of the applied light intensity (W/cm?) which range between 0 to 100 W/cm?. The responsivity 
Rypr(A/W) of the optimized structure A is around 7 A/W for A=1310 nm, and it is around 3.8 A/W for A=1550 nm. 
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Figure 16. Collector current as a function of light 
intensity at A=1310 nm and V,,=1.6 V 


B. Responsivity of the optimized structure B 
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Figure 17. Collector current as a function of light 
intensity at A=1550 nm and V,,=1.6 V 


As shown in Figures 18 and 19, the collector current I; varies as a function of the light intensity in 
phototransistor mode at Y%~=1.6 V respectively for A=1310 nm and A=1550 nm. We applied light intensity 
(W/cm?) ranging from 0 to 100 W/cm?. We noticed that the collector current I, is increasing with the increase 
of the light intensity (W/cm?). The responsivity Rypr(A/W) of the optimized structure B is around 13.1 A/W 
for A=1310 nm, and it is around 10.9 A/W for A=1550 nm. 
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Figure 18. Collector current as a function of light 
intensity at A=1310 nm and V,,=1.6 V 


C. Responsivity of the optimized structure C 
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Figure 19. Collector current as a function of light 
intensity at A=1550 nm and V,,=1.6 V 


From Figures 20 and 21, the collector current varies according to the light intensity in phototransistor 
mode at V,.=1.6 V respectively for A=1310 nm and A=1550 nm. The collector current increases gradually with 
increasing application of light intensities ranging from 0 to 100 W/cm?. The responsivity Ryp7-(A/W) of the 
optimized structure C is around 29.3 A/W for A=1310 nm, and it is around 24.3 A/W for A=1550 nm. 
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Figure 20. Collector current as a function of light 
intensity at A=1310 nm and V,,=1.6 V 


40 60 80 
Light intensity (W/cm?) 


100 


0.0016 —— 


0.0015 


0.0014 


0.0013 — 


0.0012 


Collector current (A) 


0.0011 1} -____} —__} -___} —__} ___} 
0 20 40 60 80 100 
Light intensity (W/cm?) 


Figure 21. Collector current as a function of light 
intensity at A=1550 nm and V,=1.6 V 
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D. Responsivity of the optimized structure D 

According to Figures 22 and 23, the collector current evolves in a linear way with the application of 
the light intensity in phototransistor mode at V,.=1.6 V respectively for A=1310 nm and A=1550 nm. The light 
intensity varies from 0 to 100 W/cm?. The responsivity Ryp7(A/W) of the optimized structure D is equal to 
57.8 A/W for A=1310 nm, and it is of the order of 48.1 A/W for A=1550 nm. 
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3.1.4. Comparison with other research works 

Table 3 presents the comparison between our results and some other research works, we were able to 
design photodetector topologies based on HPT InP/InGaAs. We obtained promising results in terms of 
responsivity which is a performance parameter allowing the characterization of photodetectors. These 
optimized structures can be used for digital fiber and microwave telecommunications applications. 


Table 3. Comparison with other research results on InP/InGaAs HPTs 


Reference Rupr Q=1L31 um) Rypr A=1.55 um) 
Jeong et al. [24] - 5 
Our work before optimization 3.2 2.7 
Optimized structure A 7 3.8 
Optimized structure B 13.1 10.9 
Optimized structure C 29.3 24.3 
Optimized structure D 57.8 48.1 
Ouchrif et al. [25] 10 8.4 
Khan et al. [5] 11 8.7 
HPT without optical waveguide [26] - 4 
HPT with integrated optical waveguide [26] - 29 


4. CONCLUSION 

In this present work, we investigated the impact of technological parameters on responsivity of 
InP/InGaAs HPT. Our results are mainly based on this investigation, and our objective is to design optimized 
structures using technology computer aided-design (TCAD)-Silvaco at the two wavelengths A=1310 nm and 
A4=1550 nm. We proposed four structures. The optimized structure A has a responsivity of 7 A/W for A=1310 
nm and 3.8 A/W for A=1550 nm. The optimized structure B has a responsivity of 13.1 A/W for 4=1310 nm and 
10.9 A/W for A=1550 nm. The optimized structure C has a responsivity of 29.3 A/W for A=1310 nm and 24.3 
A/W for 4=1550 nm. The optimized structure D has a responsivity of 57.8 A/W for A=1310 nm and 48.1 A/W 
for A=1550 nm. The proposed phototransistors were compared with other research works. These optimized 
structures can be integrated as photodetectors in different optoelectronic applications. 
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